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Abstract 

For the universality class of three-dimensional Ising systems the ratio of the high- 
and low-temperature amplitudes for the correlation length and for the susceptibil- 
ity are universal quantities. They can be calculated by renormalized perturbation 
theory for scalar </> theory in fixed dimensions D = 3 in the symmetric phase and 
in the phase of broken symmetry. In this article the amplitude ratios are calculated 
in the three-loop approximation. Using the fixed point values of the coupling con- 
stants we obtain /+//_ = 2.013(28) and C+/C_ = 4.72(17). 

PACS numbers: 05.70.Jk, 64.60.Fr, ll.10.Kk, 05.50.+q 
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1 Introduction 

Field-theoretic methods are a standard tool to obtain quantitative results about the 
physics of second-order phase transitions. In particular they have been applied to the 
determination of universal quantities like critical exponents or ratios of critical ampli- 
tudes. These take the same values within large universality classes, characterized by the 
dimensionality D of space and the number n of components of the order parameter. 

^Present address: Institut fiir Theoretische Physik E, RWTH Aachen, D-52056 Aachen 
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A quantity of theoretical and experimental interest is the ratio of correlation length 
amplitudes /+//_. It is defined by the behaviour of the correlation length £ as a function 
of the temperature T near the critical temperature T c through 



, f+t~ v , t > T-T c 
£~ ( J+ ' ; t := -. 1 



In this article we consider this amplitude ratio (as well as the susceptibility amplitude 
ratio) for the universality class of the three-dimensional Ising model. The corresponding 
field-theoretical model, which is assumed to be in the same class, is 4 theory in D = 3 
with a one-component (n = 1) real scalar field. The amplitude ratio can be calculated 
by means of renormalized perturbation theory in fixed dimensions D = 3. In a previous 
work M the calculation was done up to two loops. For the theoretical background and 
more details of the method we refer to that article and the references cited therein. In 
the present work we have extended the calculation to third order of perturbation theory 
and introduced a new scheme to reduce the relevant number of Feynman graphs. 
The lagrangian density for the symmetric phase (t > 0) is given by 

C(M = i(d0 o +(aO) 2 + V(0 o +) 

VOo+) = -m 2 0+ (f)l + (x) + -g $ + (x). (2) 
For the phase of broken symmetry (t < 0) we use a shift in the field variable and 



expand the potential around the minimum value Vq := J3rriQ-/go omitting constant 
terms. 

cj) -(x) := (f)(x) - v , (3) 
£(0 O -) = ^ o -(x)) 2 + V(0o-) 

V(0 O -) = ^«io- 00-0*0 + j ] ^%m Q -4>l^{x) + ^ g <Pt-(x) ■ (4) 

Besides the 4 self-interaction there is a cubic 3 interaction. This term gives rise to a 
large number of additional Feynman graphs, many of them containing tadpole subgraphs. 
To simplify the calculation and to reduce the number of graphs, we have used a method 
to eliminate tadpole subgraphs by means of a modification of the mass and coupling 
parameters [@. 

On the three-loop level no new divergencies arise. So we only have to handle subdi- 
vergencies, which are isolated by dimensional regularization. 

In the next section we summarize the renormalization scheme. The results of the 
perturbative calculations are presented in section 3, and the numerical estimates are 
discussed in section 4. 
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2 Renormalized perturbation theory 



The correlation length considered here is defined through the second moment of the 
connected two-point function. 

1 Jd D xx 2 G^(x) 



dp 



2D jd D xG. 



(2,0), 



X 



<&%) 



p=0 



Gf fi \x) := (0o(x)0o(O)) - (0o(x))(0o(O)> . 
The two-point function is related to the two-point vertex function I^ 2 ' ) by 



-r, 



(2,0), 



( P ) = (G^(p) 



-l 



(5) 
(6) 

(7) 



The renormalized mass rriR is defined by 



m 



lff 0) (p) 



p=0 



1 

e 



(8) 



and coincides with the inverse of the correlation length. 



As -r 



(2,0) 



m o + P 2 ~ where is the sum of all one-particle irreducible two- 
point graphs with amputated external legs, we can determine the correlation length £ 
perturbatively. 

The renormalization constants are defined through 



Zs 1 :- 



drr\p;m , 90 ) 



dp 2 



and 



-^({OjO};^,^) 



d 



p=0 



rj, 2 ' 0) (0;m ^ ). 



drriQ 



(9) 



(10) 



In order to obtain the universal amplitude ratio /+ / /_ various quantities have to be 
expanded in powers of a dimensionless renormalized coupling ur. The following renor- 
malization schemes are used in the two different phases. 

Symmetric phase 

The renormalized coupling constant in the symmetric phase g$ is, following |3|, [|, defined 
by the value of the four-point vertex function at zero external momenta. 



(4) 

9r ■-- 



z!r^°\{oy, mo ,g Q ). 



(11) 



It is related to the dimensionless renormalized coupling ur by 



(4) 4-D 
9R = m R U R 



(12) 
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Additionally we define the coupling renormalization constant 



-1 



lf' 0) ({0};m ,So) 

9o 



The complete set of renormalization conditions are 

r§' 0) (0; m R , u R ) = 

p(2,0 

dp 2 



-m 



R 



-1 



p*=0 

r { i' 0) ({0};m R ,u R ) 



(4) 

-9r 
-1. 



(13) 

(14a) 
(14b) 

(14c) 
(14d) 



Phase of broken symmetry 

In the phase of broken symmetry the field has a non-vanishing expectation value 



v = v + G { c' 0) , where v = y3m§/</ . (15) 
The renormalized field expectation value is 

v R := -j=v. (16) 



We define the renormalized coupling constant gR through the one-point function according 
to § 

9R--=~f- (17) 
The dimensionless coupling ur is again introduced by 

9R = m R ~ D u R . (18) 

The renormalization scheme in this phase is thus summarized by 

r2' 0) (0; m R ,u R ) = -m\ (19a) 

r (2,0 

dp 2 



9 -iffe rn R ,u R ) 



p 2 =0 



(19b) 



3m 2 

= mjf^ = g R (19c) 



rf^O};™*,^) = -1. (19d) 
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3 Perturbation series 



As the method was already developed in |jj the main problem in the present calculation 
was the number of diagrams in third order of perturbation theory. According to the renor- 
malization conditions ( [L4a| -d, |19a[ -d) we had to calculate the two-point vertex function 
and its momentum derivative at zero momentum in both phases. In the symmetric phase 
the four-point function at zero momenta and in the phase of broken symmetry the vacuum 
expectation value of the field had to be calculated, too. From these series the masses and 
coupling constants in both phases are derived. In order to distinguish the parameters in 
the two phases we label them with an index + for the symmetric (high temperature) and 
— for the broken-symmetric (low temperature) phase. The natural expansion variables 
are the dimensionless renormalized couplings 



(4) 

9r± 

,4-D 
R+ 



m 



and Ur. 



9r- 



m, 



,4-D 



The bare dimensionless coupling is defined analogously with the bare parameters 

9o 



MO 



rn, 



4-D 



(20) 



(21) 



The number of diagrams we encountered is already non-negligible: there are 204 
one-particle irreducible diagrams contributing to the inverse propagator at the three- 
loop level, compared to 20 at two loops. After exploiting symmetries there are still 162 of 
them. Many of them contain tadpoles and their elimination by means of Dyson- Schwinger 
equations simplifies the book-keeping very much ||. The reduced set of one-particle irre- 
ducible propagator-diagrams without tadpoles only contains 34 elements. Nevertheless we 
checked the calculation by means of the usual perturbation theory. The program QGRAF 
by P. Nogueira was helpful in verifying the completeness of our list of diagrams. 

The starting point of the calculation are the expansions of the renormalized masses 
and couplings in terms of the bare coupling. For this purpose a regularization scheme has 
to be used. The final results are independent of the choice of the regularization scheme. 
We decided to use dimensional regularization with D = 3 — e. 

In the symmetric phase the expansions are 



m 



R+ 



(4) 

9r+ 



m 

+ 

9o 
+ 



1 " IT i 1 ~ § (7 + ln V - 2 ) + °(t 2 )) +(^r~ ~B dlV 



8vr 

131 71. 4 32 l ndAe 

1 In - H a + -BT V 

216 27 3 3 6 1 



V162 



Mo 

8tt 



8tt) 



+ o (4) 



(22) 



|£(l-i(7 + ln^) + ^) + (-^ + 2(l + 0(e )))(g 



r99__ 373 ln 4_128 Tet _l 
1296 54 3 3 4 1 



^ 3 + °« 



(23) 
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UR+ 



Uq 



+ 



l-^(l-f(7 + ln^ + 2) + 0(e 2 )) + 



/329 1 
V216 + 6 



ZB div ] ( ^ 



u 



in 



9 9 3 



div 



8tt 



(24) 



To first order in the coupling constant Uq we have to keep terms of order e, because in sec- 
ond order there are divergent terms B dm which are proportional to e~ l . By multiplication 
they give finite contributions in the third order. 

The divergent terms B dtv and Bf iv cancel out in the final results and need not be 
displayed here. 7 is Euler's constant, and the other constants used in these equations are: 

a = g - I In 2 (I) - 1 In (§) - l L i 2 (i) = 0.0324645 (25) 



with the dilogarithm 



U 2 (x) = - 



* , In 1 - t) 

dt— - 

t 



(26) 



and 



C Tet = 7T- 6 J d z k l d i k 2 d i hA(k l )A(k 2 )A(k 3 )A{k l - k 2 )A(k 2 - k 3 )A(k 3 - h) 

= 0.1739006, (27) 

where A(k) = (k 2 + 1) . This is the only integral which we could not solve analyti- 
cally. The numerical value stems from [Bl graph 12U4] and was confirmed by our own 
calculations. 

In the broken-symmetry phase we get 

3 u 



m 



R- 



9r- 



UR- 



m 
+ 

go 
+ 

Uq 
+ 



1 + sS( 1 -i(^ + 1 ' 1 ^- 10 ) + o ( £2 ») + 



S7T 

101245 



/39T3 2_ Bdlv 
V5184 3 



8tt 



21535, 4 
H ln- 

41472 2592 3 



1723 
IT' 



3345 
1024' 



C 



Tet 



div 



MO V 
8tJ 



0(4) 



(28) 



1 - 



7% 
4 8vr 



(l-f( 7 + ln^-l)+(9( e 2 



+ 



17099 (Uq 



4051 21319, 4 1045 2849 

1 In a 

576 1296 3 24 512 



5184 V8tt 



(29) 



l_51l* (i_§ U + l n ^^) +0 (e' 

irstt v 2V' 1 7T 31 y v 



16 87T 

284719 21247, 4 

1 In - 

73728 1728 3 



( 40957 _ l B div 

V 13824 3 



8tJ 



819 

~32 a ~ 2048 



8051 Tei 31 ^ 
24 1 



Mo Y 
8n) 



(30) 



For the calculation of the universal amplitude ratio of correlation lengths we need the 
functions 



F±(u R ±) :-- 



dm 2 R± 
dm 2 ± 



(31) 



.mi 



When the bare masses and couplings are expressed in terms of the renormalized ones the 
perturbation series are 



f + (u r+ ) = i-i^±_ir^±v 

+ y H+J 2 8tt 6 V 8tt J 



+ (~^-^)(^) (32) 



27 54 3 3 J V 8tt 
3 233 / 



16 8tt 768 V 8tt 
297256 21535 

663552 ~ 5184 r " 3 ' 9(5 2048 " 7 V 8vr 



(33) 



At this point we still have two coupling constants, ur + and ur_. Following |0J we 
introduce a new coupling wr, which is defined in both phases, such that the correspond- 
ing /3-functions in both phases and consequently the numerical values of the fixed point 
couplings are equal. Previous experience suggests to choose ur such that it coincides 
with the usual coupling in the low-temperature phase: ur = ur-. Therefore we generally 
denote it by ur_ in the following. The other coupling constant expressed as a series in 
Ur- is 



Ur + {Ur_) = Ur 



1 Mr_ 1633 f Ur. 
1 + ' 



4 8vr 2592 V 8tt 

/10U239 _ 30271 4 7 2337 s f^)\ {< 
\ 165888 1296 3 8 512 J \ 8tt J V R ' 



(34) 



This relation allows to express all renormalized perturbation series in terms of a single 
coupling constant In particular for the ratio F_/F + we obtain 

F-iuRr) 



<$>-{UR- 



F + {u R+ ) 

. , 11 Mr- , 85 fu R . 



16 8vr 256 V 8vr 
/_ 109217 _ 14719 4 745 3346 N /u^y , } 



(35) 



2<&_(«U , (36) 



This function finally yields the desired amplitude ratio via 

!± 

f- 

where u* R _ is the fixed point value of the coupling and v is the correlation length exponent. 
Both u* R _ and v can be obtained in perturbation theory too, but more precise values are 
available in the literature and we shall make use of them. 
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4 Numerical results 



What is needed for the amplitude ratio of correlation lengths is the value of at 
ur- = u* R _. For the fixed point u* R _ we take an estimate u* R _ = 14.73(14) from low- 
temperature series || and another estimate u* R _ = 15.1(1.3) used in |IIJ. For comparison 
we note that the zero of the 3-loop /3-function is located at u* R _ = 14.2, using a (2,1)-Pade- 
Borel approximation. The entries in table 1 result from the possible Pade approximants 



Fixed point 


*-(«*-) 


[3,0]-Pade [2,1]-Pade [1,2]-Pade [0,3]-Pade 


14.73 | 


1 


1.5288 1.5301 1.5002 1.5149 


15.1 


10 




1.5456 1.5471 1.5133 1.5301 



Table 1: $_ as a function of the low-temperature fixed point u* E 



to $_(%_) evaluated at u* R _. The values are quite close together and the mean value is 
1.526(26), where the error represents the maximal deviation. The numerical convergence 
of the series at the fixed point is rather good, 

$_(«/?-) = 1 + 0.410 + 0.118 + 0.012 , (37) 

although it is expected to be asymptotic only. An application of the usual Pade-Borel 
summation method does not improve the result. On the contrary the uncertainty is more 
than doubled. 

Using the high-temperature coupling (u* R+ ~ 24) as an expansion parameter is much 
worse. For $ + (u/? + ) we get a mean value of 1.39(34) with Pade, and 1.48(24) with Pade- 
Borel summation, respectively, so the maximal deviation is more than ten times higher 
as above. This is related to the poor numerical convergence: 

$+(u* R+ ) = 1 + 0.657 + 0.146 - 0.396 . (38) 

Therefore we consider the estimate from $_ as more reliable. 

For the critical exponent v we used Monte Carlo results [y = 0.624(2)) |TT] and values 



of renormalized perturbation theory [y = 0.6300) [[12]. After exponentiation with these 
values for v we get the universal amplitude ratio for the correlation length £ according to 

f 

j- = 2.013(28) . (39) 

Using the high-temperature coupling instead we would get a value of 1.98(20) employ- 
ing Pade-Borel approximations. The error is ten times larger than in (j39|). As discussed 
above the low-temperature coupling appears to be the better expansion parameter. 
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With a similar method we have calculated the universal amplitude ratio of the sus- 
ceptibility. As before we get best results with the low-temperature coupling, namely 

-± = 4.72(17) . (40) 

The same ratio has been calculated by means of three-dimensional perturbation theory 
in JTJ with the result C + /C_ = 4.77(30). 



5 Conclusion 



Our third order calculation of the amplitude ratio of correlation lengths (^) is a confir- 
mation and improvement of the two-loop result of 2.03(4) Q. Theoretical estimates in 
the e-expansion (1.91 (14j]), high- and low-temperature expansions (1.96(1) [|llj, 1.94(3) 
0) are lower. Experimental values (2.05(22), 2.22(5) @, 1.9(2), 2.0(4) |7|) and Monte 
Carlo results (2.06(1) (T^|) are above or close to our results. 
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